Constraints on a variable dark energy model with recent 

observations 



Puxun Wu a ' and Hongwei Yu a '* 
a Department of Physics and Institute of Physics, 
Hunan Normal University, Changsha, Hunan 4 10081, China 
h School of Sciences and Institute of Physics, 
Central South University of Forestry and Technology, Changsha, Hunan 4 1 0004, China 

Abstract 

We place, by the maximum likelihood method, constraints on a variable dark energy model with 
the equation of state w = wq/[1 + 61n(l + z)] 2 using some recent observational data, including 
the new Sne la data from the SNLS, the size of baryonic acoustic oscillation peak from SDSS 
and the CMB data from WMAP3. We find that the SNLS data favor models with wq around 
— 1, in contrast to the Gold data set which favors a more negative wq. By combining these three 
databases, we obtain that tt m = 0.27tomt, ™o = ~ 1 - ll -oio and b = °- 31 ia3i with X 2 = H0.4 at 
the 95% confidence level. Our result suggests that a varying dark energy model and a crossing of 
the w = — 1 line are favored, and the present value of the equation of state of dark energy is very 
likely less than —1. 

PACS numbers: 98.80.-k, 98.80.Es 
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is not dnven_b y dark energy such as modffied gravrty H , theories with compacted extra 

9] and Cardassian cosmology [10( etc. 



I. INTRODUCTION 

Since the Hubble diagram of Type la Supernovae (Sne la) [lj first indicated that the 
universe is undergoing an accelerating expansion, many works have been done, trying to 
explain this phenomenon. A large number of models for the cosmic acceleration are proposed 
by assuming the existence of an energy component with negative pressure in the universe, 
named dark energy, which at late times dominates the total energy density of the universe 
and drives its acceleration of expansion. The simplest candidate of dark energy is the 
cosmological constant A Q], with the equation of state w = p/p = —1. However, the 
cosmological constant has to be extremely fine-tuned in order to induce the observed cosmic 
expansion, and this has led many authors to use scalar fields, such as quintessence j^], 
phantom [4|, quintom jlj] and Chaplygin gas as alternative models for dark energy. On 
the other hand there are also some other models where the observed cosmic acceleration 

a 

dimensions |8|], DGP model 

In order to determine the recent expansion history of our universe, one can also use 
a different approach, that is, to assume an arbitrary parametrization for the equation of 
state w(z) for dark energy, where z is the redshift. The parametrization may not be moti- 
vated by any particular fundamental physical theory and is thus "model-independent". It 
however needs to be designed to give a good fit to the observational data. The simplest 
parametrization is w = constant. Some other proposals, including w(z) = wo + w\z [llf . 
w(z) = w + W\z/{\ + z) etc, are also made. Recently Wetterich proposed an inter- 
esting phenomenological parametrization for a variable dark energy, in which the effective 
equation of state is expressed as: 

w{z) = [i + bHi + zW ' (1) 

where w represents the present value of the equation of state and b is a positive constant 
characterizing the change of w(z) with redshift. Apparently with this w(z) there are three 
model parameters (Q m ,w ,b) to be determined by the observations, where fl m = p m /p c 
represents the present matter density parameter and p c = 3Hq/8ttG is the present critical 
density of our universe. The advantage of this parametrization over other parameterizations 
of a time varying equation of state, w(z) is that it covers the whole available redshift range 



while other parameterizations proposed before cover only a restricted range of redshift 
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Later, using Cosmic Microwave Background, Large Scale Structure, and Sne la 



data, Doran, Karwan and Wetterich 



15| discussed the wq and the dark energy fraction at 



very high redshift Q d in this model, and found at the 95% confidence level Wq < —0.8 and 
Vt e d < 0.03, where Vt e d and wq are related to b by b = — 3wo(ln ^-j^ + ln -tw 1 ) -1 with Q d being 

d d 

the present energy density of dark energy. Movahed and Rahvar 16] have used the Gold 
Sna la data [isl |. the position of first acoustic peak of the Cosmic Microwave Background 
radiation (CMB) and the size of baryonic acoustic oscillations peak to constrain this model 
and obtained at the 2a confidence level fl m = 0.27+q;q3, w 



-1.45l°;f and b 



Recently Asitier et al. 19] released the data of high redshift supernovae from the Su- 



pernova Legacy Survey (SNLS). In this survey the systematic uncertainties and systematic 
errors are reduced. It is worth noting that the SNLS data set is a better agreement with the 
WMAP data compared to the Gold Sne la set 17|. Thus, in this paper we will reexamine 
this variable dark energy by using the 115 new SNLS Sne la data, the size of baryonic acous- 
tic oscillations peak detected in the lar ge-s cale correlation function of luminous red galaxies 



from Sloan Digita 



WMAP result 
and b = 0.31 



Sky Survey (SDSS) 20j and the CMB data obtained from the three-year 



2jJ. We obtain at a 95% confidence level Q m = 0.27^;^, w = -1.11 



1+0.21 
-0.30 



-0.71 
-0.3+ 



II. THE BASIC EQUATION 



By using Eq. and the equation of energy conservation, it is easy to obtain the evolution 



of density of dark energy 
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16] 



Pd = Pdo(l + z 



.3(1+0(2)) 



(2) 



where p^o denotes the present density of dark energy and w(z) 



, , , , , Since WMAP 

l+6ln(l+z) p-i 



observations strongly indicate that the geometry of our universe is spatially flat [22j, we 
will ignore the term containing curvature factor in Friedman equation. If the radiation 
components in universe are further ignored, we find for the Hubble parameter 



H 2 (z; tt m , w ) = H$[Q m (l + zf + (1 - Q m )(l + z 



,3(1+15(2))! 



(3) 
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where Ho is the present Hubble constant. Meanwhile one can show that for a flat universe 
the Luminosity distance, d L , can be expressed as 

r r z d?' 

d L (z,H ,Q m ,w ) = —(l + z) —— s . (4) 

H Jo E{z',VL m ,wo) 

Here c is the velocity of light and E(z, Q m , Wq) = H(z; fl m , w )/H . 



III. CONSTRAINTS FROM SNE IA, SDSS AND CMB DATA 



The the 115 new Sne la data includes 44 previously published nearby Sne la and 71 
distant Sne la released recently by the Supernova Legacy Survey (SNLS) 19] which is a 
planned five year survey of SNe la with z < 1. Constraints from Sne la can be obtained by 
fitting the distance modulus n(z) 

fi(z) = 5lo gl o[D L (z)]+M . (5) 

Here D L = H Q d L and M. = M — 5logio(H ) , M being the absolute magnitude of the object. 



Recently Eisenstein et al 20j successfully found the size of baryonic acoustic oscillation 
peak using a large spectroscopic sample of luminous red galaxy from the SDSS and obtained 
a parameter A, which is independent of dark energy models and for a flat universe can be 
expressed as 
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dz 



2/3 



A EizS^Ui E(z)\ ' (6) 
where Z\ = 0.35 and A is measured to be A = 0.469 ± 0.017. Using parameter A we can 
obtain the constraint on dark energy models from the SDSS. 

For the CMB data, the shift parameter R can be used to constrain the dark energy models 
and it can be expressed as [23j 

dz 



R 



E(z) 



(7) 



for a flat universe, where z r = 1089. From the three- year WMAP result 221 ]. the shift 



2i 



parameter is constrained to be R = 1.70 ± 0.03 

In order to place limits on model parameters (Q m , wo, b) with the observation data, we 
make use of the maximum likelihood method, that is, the best fit values for these parameters 
can be determined by minimizing 

,2 



X 



v \pLobs(zi) - fJ,{zi)} 2 (A- 0.469) 2 (R-1.70f 

n 1 z z . _„ r 



0.017 2 



0.03 2 



(8) 



For the SNLS Sne la data set, at a 95.4% confidence level we obtain Q m = 0.27l n 97 , 
Wo = —1.03^42 and b = 0.0 +4 ' 3 . These are different from the result obtained in Ref. [16| 
using 157 Gold Sne la data, where VL m = 0.01+^, w = -1.90t 3 H and b = 6.00^. 
Apparently the SNLS data set favors models with wq around —1 while the Gold set favors 
a more negative wq. Meanwhile the best fit value of b(= 0) for SNLS data show that a 
non-varying dark energy model is favored. However, if combining the SNLS Sne la, SDSS 
and CMB, we find that at a 95% confidence level tt m = 0.27^!, = -l.lllo.iS and 
b = 0.31 

-oil with x 2 — H0.4. The best fit values show that a variable dark energy model 
is favored since b is nonzero and it is very likely the present value of the equation of state is 
less than —1. In Fig. ([1]) the Hubble Diagram for 115 SNLS Sne la data set is shown with 
the (fi m , wo, b) = (0.27, —1.11, 0.31). The contour plots of fl m and wq by fixing b at its best 
fit value 0.31 are shown in Fig. §2§. The contour plots of fl m and b by fixing wq at its best 
fit value —1.11 are shown in Fig. ([3]). In Figs, (jlj) we give the evolutionary curves w(z) vs z 
with la error bar based on SNLS Sne Ia+SDSS+CMB data. This figure shows graphically 
that SNLS Sne la +SDSS+CMB data favor a varying dark energy model. Comparing with 
the results obtained in Ref. 1(|, we find that in our results at the 95% confidence level 
stronger constraints on wq and b are obtained. Meanwhile it is easy to see that we also 
obtained a stronger constraint on wq than that obtained in Ref. [15J . 



IV. CONCLUSION 

In this paper we have placed constraints on a parameterized dark energy model [3| using 
the new SNLS Sne la data sets, the size of the baryonic acoustic oscillation peak from SDSS 
and the shift parameter from the CMB observation. It is found that a non-varying phantom 
dark energy model is favored if the SNLS Sne la data set is used in contrast to a varying 
dark energy when the Gold Sne la data is utilized 16] ; and different from the case of the 
Gold set, the model with wo around —1 is favored by SNLS data at a 95.4% confidence 
level. Combing three databases (SNLS Sne la, SDSS and CMB), we obtain a constraint on 
the model parameters (Q m ,wo, b) , which suggests that a varying dark energy model and a 



crossing of the w — — 1 line are favored 
of dark energy is very likely less than — 1. 



24|, and the present value of the equation of state 
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FIG. 1: The Hubble diagram for 115 SNLS Sne la data with the best fit parameters (£l m ,wo,b) 
(0.27, -1.11,0.31) obtained from the combination of SNLS, SDSS and CMB databases. 
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FIG. 2: The la, 2cr and 3a confidence contours for Q m and wq with 6 at its best fit value 0.31 from 
the combination of SNLS, SDSS and CMB databases. 



9 



0.8 



0.6 



x> 0.4 



0.2 




w ( j=- 1.11 



/ 1 

0.1 0.2 0.3 0.4 0.5 



FIG. 3: The la, 2a and 3d confidence contours for £l m and b with wo at its best fit value —1.11 
from the combination of Gold, SNLS, SDSS and CMB databases. 
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FIG. 4: The behavior of w(z). The solid line plots w(z) by using the best fit parameters (wq = 
-1.11,6 = 0.31) obtained from the SNLS+SDSS+CMB data and the dotted lines are for la errors. 
The solid line shows clearly a varying equation of state parameter. 
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